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! ABSTRACT 

This is the second of two papers investigating the spherical averaging of eUipsoidal 
galaxy clusters in the context of X-ray and Sunyaev-Zel'dovich (SZ) observations. In 
the present study we quantify the orientation-average bias and scatter in observables 
that result from spherically averaging clusters described by ellipsoidal generalizations 
of the NFW profile or a nearly scale-free logarithmic potential. Although the mean bi- 
(~| ] ases are small and mostly < 1%, the flattest cluster models generally have a significant 

C^' mean bias; i.e., averaging over all orientations does not always eliminate projection 

^ I biases. Substantial biases can result from different viewing orientations, where the 

i-^ • integrated Compton-y parameter (Fsz) and the concentration have the largest scat- 

ty I ter (as large as tr ~ 10% for Fsz), and the emission-weighted temperature (Tx) has 

the smallest (ct<0.5%). The very small scatter for Tx leads to Yx and Mgas having 
virtually the same orientation biases. Substantial scatter is expected for individual 
CsJ ■ clusters (up to a 8%) in the correlation between Ysz and Yx in comparison to the 

\ small mean bias (a < 1%) applicable to a random sample of clusters of sufficient size. 

■ For ellipsoidal NFW models we show that the orientation bias for the total cluster 

, mass attains a minimum near the radius r25oo so that the spherically averaged mass 

\^ • computed at this radius is always within « 0.5% of the true value for any orientation. 

Finally, to facilitate the accounting for orientation bias in X-ray and SZ cluster stud- 
' ies, we provide cubic polynomial approximations to the mean orientation bias and la 

scatter for each cluster observable as a function of axial ratio for the ellipsoidal NFW 
models. 
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1 INTRODUCTION 



Galaxy clusters are potent tools for cosmological studies (e.g., see reviews by iHenrvlbooal : ISchueckejlioOsI : lArnaudlliooi : IVoitI 



20051 : iTozzil 120071 lAllen et al.ll2011f ). The redshift evolution of cluster baryon fractions supplements well-known geometrical 



cosmological constraints obtained by the cosmic microwave background radiation, supernova distances, and baryon acoustic 
oscillations. Important cosmological information complementary to these probes is provided by the growth of cosmic structure, 
which is manifested in a variety of ways for clusters, notably via the evolution of the cluster mass function. Since the hot 
intracluster medium (ICM) dominates the cluster baryons. X-ray observations are essential for measuring baryon fractions. 
They also provide one of the most effective means to measure the total gravitating mass for a cluster that obeys approximate 
hydrostatic equilibrium. 

To obtain the most accurate constraints on the baryon fraction and total gravitating mass with X-ray observations 
requires obtaining accurate measurements of both the radial ICM density and temperature profiles out to a large fraction of 
the virial radius. These demanding requirements impose rather strict limits on the number of clusters which can be usefully 
analyzed with current data and data anticipated in the near future. Consequently, for cosmological studies with X-rays there is 
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much interest in using global scaling relations of quantities that act as mass proxies; e.g., the ICM temperature and gas mass 
measured from X-ray observations and the ICM pressure measured from either X-ray or Sunyaev-Zel'dovich (SZ) observations. 
Useful measurements of such proxies are possible for data of much lower quality and therefore may be applicable to much 
larger cluster samples than if detailed radial ICM profiles are needed. Moreover, even if the measurements of such mass proxies 
are of low quality, by averaging over large samples the impact of large statistical errors on measurements for individual systems 
is mitigated leading to precise cosmological constraints, which are accurate provided systematic errors also "average out" or 
are understood and controlled. 

The hot ICM also traces feedback processes during cluster for mation and evolution. For purely gravita tional evolution 



the radial ICM entropy profile behaves as approximately r ' (e.g., iTozzi fc NormanI l200ll : IVoit et ah 2005 ), whereas feed 



back from an AGN and supernovas flattens the entropy profile mostly in the central regions (e.g., Tozzi fc NormanI 



2008 



- , gp -- ., - „ 

Brighenti fc Mathewsll200ll : ]voit et al. 200 2: Bor gani et alboOSl : iRomeo et ahlbood : IVounger fc BrvanlboO?! : [McCarthy et al 



2001 



2010|). Seemingly at odds with these theoretical expectatio ns, recent X-ray studie s of the Perseus an d Virgo clus- 



ters measured en t ropy profiles that flatten near the virial radius (jSimionescu et al.ll201ll : lUrban et al.l 120111 ). In addition. 



Simionescu et al.l l|2011f ) measure a gas fraction for Perseus that rises well above (over 50%) the cosmic mean baryon frac- 



tion. To reconcile these unanticipated results with the theoretical picture mentioned above, the authors of these studies 
postulate that near the virial radius the cluster ICM is substantially multiphase due to clumping from incompletely mixed 
materi al currently falling into the cluster - a hypothesis that subsequently has been given support by a cosmological simu- 
lation (jNagai fc LauluOllI ). In light of these results it is interesting that our recent Chandra and Suzaku study of the fossil 
group/cluster RXJ 1159-1-5531 out to the virial radius does not find a flat entropy profile an d also yields a measurement of the 
baryon fraction that is very consistent with the cosmic mean value IjHumphrev et al.ll2011f ). A possible interpretation of these 
different results is that the magnitude of the gas clumping varies between clusters such that the most relaxed systems (like 
fossil groups) display the least amount of clumping. Of course, verifying this interpretation requires an accurate accounting 
of the important systematic errors involved with such measurements. 

Deviations from hydrostatic equilibrium likely account for the largest systematic errors in X-ray measurements of cluster 
masses. Theoretical studies of clusters formed in hydrodynamical cosmological simulations find that such deviations can lead 
to mass errors ranging from a few percent up to ^ 30% even for nominally relaxed clusters (e.g., Tsai et al]|l994 : Evrard et al 



19961 : iRasia et allbood : iNagai et ahlboOTl : IPiffaretti fc Valdarninilboosi : I Vazza et allboill ). The simulations typically find that 



non-thermal pressure generated by turbulent motions in the ICM is primarily responsible for the systematic underestimates 
of the masses. Whi le the modest turbulent velocities associate d with these mass biases are not strongly const rained by current 
spectral data (e.g., Schuecker et al. 20041 : Sanders et al. 2011), the calorimeter to be operated on ASTRO-H (|Takahashill201ll ) 
should measure the amount of turbulent motions with novel precision, thereby providing a direct constraint on the accuracy 
of the hydrostatic equilibrium approximation and substantially reducing this source of systematic error for many clusters. 

Arguably next in importance is the systematic error associated with the ubiquitous assumption of spherical symmetry, 
since it is well-known that clusters are not spherical. While this assumption is not required, it is the default standard for 
several reasons. First, the intrinsic shape and orientation of a cluster, even if it is assumed to be a triaxial ellipsoid, are 
not directly observable. Second, although promisi ng approaches exist to constrain these properties indirectly though various 
combinations of X-ray, SZ, and lensing data (e.g., Zaroubi et al. 19981 ). they necessarily require more data and involve more 
complex modeling than if spherical symmetry is assumed. Finally, while it is generally recognized that for an individual 
cluster substantial error might be expected by assuming spherical symmetry (e.g., when viewing a highly flattened spheroidal 
cluster face-on), it is widely assumed (or hoped) that if it were possible to average such measurements over all possible 
cluster orientations, or rather average results over a statistically large cluster sample, that the effects of random projection 
orientations on cluster measurements should be eliminated entirely. 

Over the years several theoretical studies have mentioned that observed IC M proper ties should not va ry dramatically as 
a result of different intrinsic shapes and orientations. The pioneering study by Binnev fc Strimpe| j 19781 ) showed that even 
for spheroidal clusters highly inclined to the line-of-sight, deprojection assuming that the symmetry axis lies in the plane of 
the sky (effectively an assumption of spherical symmetry for large inclination angles) did not produce very large errors in the 
deprojected X-ray emissivity proflle. It has been known since the first hydrodynamical cosmological simulation of a cluster that 
also simulated X-ray observations to quantify errors in the recovered mass, that significant 25%) errors c ould be expected 
from a combination of deviations from hydrostatic equilibrium and spherical symmetry (jTsai et al.M 19941 ). More recently. 



cluster simulations have shown that significant scatter in the measurement of the integrated Compton-y par ameter (Ysz 
(proportional to the integrated ICM press ure) arises from effects due to non-spherical cluster morphology (e.g.. IWhite et al 



(proportional to tne mtegrated iUiVi press ure) arises irom ettects due to non-spnerical cluster morpnology le.g.. IWmte et al. 
20o3 : IShaw et aLlliooi : iKrause et ahlboill) . But precisely how much scatter arises from morphology alone is difficult to extract 



from such simulations; e.g., the large s catter in the relati onship between measures of cluster morphology and deviations in 
the correlation between Ysz and mass (jKrause et al.ll201l[ ). 

It is therefore worth examining comparatively simple non-spherical models that effectively isolate the impact of assuming 
spherical symmetry on cluster measurements in X-rays or the SZ effect. Such examinations are quite uncommon in the 
literature, and when they do occur, they most often take the form of ancillary systematic error checks on X-ray cluster 
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measurements of the total mass and gas mass. For e xample, to examine the error from assuming spherical symmetry on their 
measurements of the gas mass in the cluster A 478, IWhite et al.l (|l994r ) divided up the radial surface brightness profile from 
ROSAT into different azimuthal sectors, concluding that the masses were consistent between the different sectors within the 
estimated 3a statistic al errors. This approach has also been applied to other systems, including the more recent study by 
Churazov et al. I (|2008l ) of Chandra observations of the low-mass clusters, M 87 and NGC 1399, in which they conclude errors 



of a fe w pe rcent in their total mass me asurements can be attributed to the assumption of spherical symmetry. iFabricant et al. 



(1984) and Buote fc Canizare^ (1996b) noted in their studies of the intrinsic shapes of some clusters with Einstein and ROSAT 
that the spherically averaged total mass profiles obtained using both oblate and prolate sp heroidal models were consi stent, 
within the statistical errors, with the results obtained from assuming spherical symmetry. Buote fc Canizare3 ( 1996lj ) also 
mentioned small (~ 2%), but significant, errors on the gas mass if spherical symmetry is assumed. Finally, using spheroidal 
models , Grego et al.. (2000.) estimate ~ 20% systematic errors from assuming spherical symmetry on their SZ-derived estimates 
of t he gas mass in A 370. 



Piffaretti et al. I (|2003t ) and iGavazzil (|2005l ) conducted the first studies dedicated to quantifying biases owing to the as- 



sumption of sp h erical symmetry in measurements of total cluster masses, in particular with respect to X-ray observations 



Piffaretti et al. (|2003l ) presented results for ten clusters observed with ROSAT. They modeled each cluster ICM as a tri- 



axial ellip soid projected down its intermediate principal axis, with the major-axis profile described by the isothermal j3 
model (jCavaliere fc Fusco-Femiandll976 ; Sarazin,1986, ). Their models indicated that the biases from assuming spherical sym- 
metry were ^4% for the total mass and 5, 5% for the gas fraction . Gavazzi ( 20051 ) conducted a predominantly theoretical 
study of clusters modeled as spheroidal NFW ([Navarro et al.l 119971 ) mass distributions and isothermal ICM . To emphasize 
projec tion effects he focused on the special case of a cluster viewed down its symmetry axis (i.e., face-on). Like lPiffaretti et al. 
( 2003h . generally he found mass biases of only a few percent owing to the assumption of spherical symmetry. 

While significant biases from assuming spherical symmetry have been clearly demonstrated for particular intrinsic shapes 
and viewing orientations, the critical question of whether such biases completely "average out" when viewed from all directions 
remains unanswered. To address this important issue we compute theoretical orientation-average biases and scatter for ideal 
X-ray and SZ measurements as a function of the intrinsic shape of a cluster. We also consider several observables of interest 
to X-ray and SZ studies - concentration, total mass, gas mass, gas fraction, emission- weighted temperature, Ix, and Isz - 
and a wider range of models, including models with ICM temperature gradients. Finally, for a subset of models we provide 
cubic polynomial approximations to the mean orientation bias and la scatter of each observable that may be convenient for 
assessing the errors from assuming spherical symmetry in X-ray and SZ cluster studies 

The paper is organized as follows. We define the ellipsoidal models in [J2] and specify the default parameter values in [J3l 
In [)4]we describe the procedures to produce theoretical and observed spherically averaged observables. We present the results 
in !j5]and our conclusions in [JG] Our computations employ several of the analytical results for the deprojection and spherical 
averaging of ellipsoidal potentials derived in our companion paper (jBuote fc Humphrevll201ll . hereafter Paper 1). 



2 MASS MODELS 



It is customary to approximate the isodensit y contours of cluster dark matter halos formed in cosmological simulations by ellip- 
soids of constant shape and orientation (e.g., Jing fc Sutol20o3 : Bailin fc Steinmet J2005 : AUgood et ahlfioO^ : Munoz-Cuartas et al, 
201lh . While this is often a satisfactory approximatio n, more ge nerally the isodens ity contours of d ark matter halos exhibit 
variations of shape and orientation with radius (e.g., Bailin fc Steinmetz 2005: AUgood et al. 20061 ) . To explore a variation 
of shape with radius in the mass distribution we also cons ider models where the isop otential contours have constant shape 
and orientation; e.g . , the cored logarithmic potential (e.g., Binnev 1981 : Evan j T993I ) . Here we extend the nomenclature of 
Kassiola fc Kovneil (l993l ) to three-dimensional systems by referring to models having ellipsoidal isodensity surfaces of con- 
stant shape and orientation as "Ellipsoidal Mass Distributions" and models with ellipsoidal isopotential surfaces of constant 
shape and orientation as "Ellipsoidal Potentials", the latter of which are the focus of Paper 1. 

Since each of these model types determines only how the fiattening of th e mass distribution v aries with radius, the radial 
profile must also be specified. We will consider two different models: NFW ([Navarro et al.lll997l ) and the cored logarithmic 
potential mentioned above. The NFW profile provides a good overall description of dark matter halos formed in cosmological 
simulations an d also of the total mass profile inferred f rom X-ray observations of massive elliptical galaxies and clusters (e.g., 
see reviews bv lArnaudlliooil : iBuote fc Humphrevlboi^. and references therein). Althou gh the centers of dark matter halos are 
better represented by a Sersic/Einasto model ( Navarro et al. 2004; Merritt et al. 20051 ) . the additional accuracy at small radii 
is unimportant for our present study which focuses on global halo structure. Unlike the NFW model, the cored logarithmic 
potenti al is not motivated by cosmological simulation s. Nevertheless, it also corresponds to the widely used isothermal /3- 
model ( Cavaliere fc Fusco-Femiano 19761 : Sarazin 19861 ) which is well-known to provide a good description of the radial X-ray 
surface brightness profiles of many galaxies and clusters. 

Finally, we take the ICM to be a tracer population within a cluster, which is a reasonable approximation given that the 
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observed cluster gas fraction is ~ 10% within the radius (rsoo) we investigate (e.g., Pratt et al. 201Cll V Although we ignore 
the self-gravity of the ICM we do examine how the gas mass (and gas fraction) implied by hydrostatic equilibrium in the 
single-component potential are affected by spherical averaging (see Q. 



2.1 Spherical Models 

We begin by presenting the density, enclosed mass, and gravitational potential for a spherically symmetric system for both 
the NFW and cored logarithmic potential models. This provides an opportunity to define some useful quantities that will be 
used throughout the paper and will serve as a convenient point of comparison for the ellipsoidal generalizations defined below 
in §[23]and §[Ml 



2.1.1 NFW 

The NFW mass density profile, p(r) oc r^^{rs +r)^^, has two free parameters: (1) a normalization and (2) a scale radius (r^) 
denoting where the logarithmic slope of the density profile equals -2. It is conventional to replace these with the following 
related parameters. Define the concentration parameter, c = VA/rs, where va is defined so that the average density of the 
halo with that radius equals A times the critical density of the universe (pcrit). Typically quoted values for A range from 
100-2500. For the other parameter, define Ma to be the mass enclosed within ta. With these definitions, the NFW density, 
enclosed mass, and gravitational potential take the form. 



p{f) = 



Ma 



Anr^ A{c) cf{l + cf) 



(1) 



M(< f) = Ma 



A(c) 



ln(l -f- cf) 



1 + cf 



(2) 



$(f) = 



GMa 1 ln(l + cf) 
rA A{c) f 



(3) 



where. 



r 

TA ' 



A{c) = ln(l + c) 



1 + c 



3Ma 



Apcrit 



(4) 



2.1.2 CORELOG (also Isothermal P-model) 

As befits its name , the cored logarithmic pote ntial is defined by the form of its potential, <I> = —(110/2) ln(r^ +r'^) (e.g.. lBinnev 
198ll : lEvanj|l993l : ISinnev &: TremainelboOa ). Like the NFW model it has two free parameters: a normalization (vo) and a 
scale radius (rc), usually called a "core radius." For consistency of presentation we redefine these parameters in analogy with 
that done for the NFW model above. That is, we define a concentration parameter, c = ta/tc, and mass. Ma, where ta is 
again defined so that the average density within a sphere of that radius equals A times the critical density. When expressed 
in terms of c and Ma the cored logarithmic potential and derived quantities take the form. 



p(f) = 



^^^(l + c^) 



47rr^ 



3 {cff 
[l + (cf)^]" 



(5) 



M(<f) = Ma(1 + . 



1 + {cfY 



(6) 



<l>(f) 



GMa (1 + 
2rA 



In 



1 + 
1 + (cf ) 



(7) 



where, f and Ma, are defined as above (eqn. |4]). We have chosen the arbitrary constant term in the potential so that $ = 
when r — TA- Henceforth we will refer to this model as the "CORELOG" model. In order to provide a sharper contrast 
with the NFW model, we will focus on the case c ^ 1 so that $ ~ ln(r), corresponding to the singular isothermal sphere . 
It is noteworthy that eqn. ((7| also corresponds to the potential of the isothermal /3-model I Cavaliere fc Fusco-Femianol 19761 : 
Trinchieri et al.l 19861 ). where the ICM density is assumed to obey, pgas oc (1 -I- {r /rc)^)~^^^'' , and the factor outside the 
logarithm equals, —{3l3/2){kBT/ finia), where T is the ICM temperature, fes is Boltzmann's constant, p is the mean atomic 
weight of the ICM, and ma is the atomic mass unit. 
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2.2 Ellipsoidal Models ~ General Considerations 

As in Paper 1 we consider an ellipsoid of principal axes a, b, c with axis ratios p„ = b/a and g„ = c/a satisfying < < Pt, < 1. 
When the ellipsoid is oriented so that a lies along the x-axis, b along the y-axis, and c along the z-axis, the ellipsoidal radius 
is defined by, 

2 2 

2 2 , y Z , , 

a„=a; (8) 

Pv Hv 

For the ellipsoidal generalizations of the spherical NFW and CORELOG models we consider below, we will make use of the 
following quantities, 

. _ a„ 3Ma . 

a = — , 3" = Apcrit, (9) 

where, in analogy to the spherical case, a a is defined so that the average density within an ellipsoid a„ = a a equals Apcrit, 
and Ma refers to the mass contained within that ellipsoid. 



2.3 Ellipsoidal Potentials - EPs 



For a detailed discussion of this model type we refer the reader to Paper 1. Briefly, we define an ellipsoidal potential (EP) 
model such that the gravitational potential is an ellipsoid of constant shape q^) and orientation; i.e., $ depends only on 
a„. For this model type one begins by specifying $(a„) and then deduces the mass distribution from it. The key advantage of 
an EP is that a simple, analytic form for <1> can be adopted, enabling much faster computational evaluation than for ellipsoidal 
mass distributions ( ii2.4|) . A disadvantage is that when the potential is sufficiently fiattened (the amount depending on how 
steep is the radial potential profile), the model can be come unphysical in region s near the shortest principal axis; e.g., as 
manifested by negative values for the mass density (e.g., Binnev fc Tremainel 20081 ) . However, since we primarily make use of 
the integrated mass profile, which is as well-behaved as for spherical analogs of $ (see Paper 1), localized peculiar features in 
the density profiles of EP models are not important for our investigation. 

In Theorem 1 of Paper 1 we show that the mass enclosed within a„, A/(< a„), is related to the the potential ^{a^) in 
the same manner as the spherical M(< r) is related to $(r), with the introduction of a proportionality constant, 



v{Pv,qv) 



Pvqv 



1+ — + — 

pi qi 



(10) 



which is of order unity. We emphasize that p„ and q^ define the shape of "I> and the bounding surface for the calculation of 
M(< a„) (and thus definition of Ma below in i]2.3.1l and ^2.3.2|) . but the mass density distribution of an EP generally has 
smaller axis ratios (see 331 and is not strictly ellipsoidal. 

We now define versions of the NFW and CORELOG models appropriate for EPs. 



2.3.1 NFW-EP 



We convert the spherical NFW model to an EP starting with the potential defined by eqn. ^ and letting r ^ and 
TA — >■ flA, 

GMa 1 ln(l + ca) 



(11) 



OA A{c) a 

where c = flA/ciscaic and Oscaic is an ellipsoidal analog of r^. From this modified potential the mass profile is computed using 
eqn. (7) of Paper 1 to give. 



M(< a) = Ma 



v{pv,qv 

A{c) 



ln(l + ca) 



1 + ca 



and the density profile is computed from Poisson's equation which yields. 



p{x,y,z 



Ma 



Ana^ A{c) (ca)2 (1 + ca) 
1 + ca at \ pt qi 



pi qi 



+ 



(1 + ca) In (1 + ca) 



1 1 

1+ — + — 

pi qi 



2 1''^"'^ 4 4 

av \ Pv qv 



(12) 



(13) 



2.3.2 CORELOG-EP (also EP analog of Isothermal /3 model) 

The procedure to construct the EP analog of the CORELOG model follows the same procedure as for the NFW model above 
but starting from the potential given by eqn. ((7|). We obtain for the potential. 
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* « = — ^ — r^ln 



1 + 

1 + {caf 



(14) 



where c = aA/cic and ac is an ellipsoidal analog of rc. The mass and density profiles derived from this potential are respectively, 

~ 3 

M(< a) = 77(p„, g.)MA(l + c^) ^ " . „ , (15) 

1 + (ca) 

and, 

Ma 



(1+c^) 



1 



2 2 

2 , y ,2 

a; + — + — 



(16) 



This expression for the density is the triaxial generalization of the well-known spheroidal form (e.g., Binnev fc TremainelboO 



2.4 Ellipsoidal Mass Distributions - ElVIDs 

We define an ellipsoidal mass distribution (EMD) such that the mass density is an ellipsoid of constant shape and orientation; 
i.e., p depends only on the ellipsoidal radius a„, where now the axial ratios p„ and Qy refer to p rather than The gravitational 
poten tial generated by an EMD is not an ellipsoid but instead has the form (e.g.. IChandrasekhailll987l : iBinnev fc Tremaine 
20081 ). 



${x) — —TvGpvqvCl'A 

where, 
a [x, u) 



du 



^ ^ ^ a'^{x,u) 



da p(a 



+ 



+ 



A{u) = ^ {a\ + it)((pi,aA)^ + M)(((/„aA)^ + w). 



(17) 



(18) 



and the function A(u) should not be confused with the over-density value as si gnified by a a- Because the improper integral 
over du with upper limit oo is inconvenient for numerical evaluation, we follow iMerritt fc FridmanI l|l996l ') and define a new 
variable, 

(19) 



so that, 



a^(x, s) 



+ 



H-(ai-l)s2 1 [(p.aA)2 - l]s2 1 + [(g„aA)2 - l]s2 



A(s) = &{s) = V(l + {al - l).s2)(l + [(p„aA)2 - l]s^){l + [(g^aA)^ - l]s^ 



^(x) — — 27rGp„g„aA 



ds 



da'^ p{d''^) 



We use eqn. (|22|) for all numerical evaluations of EMD potentials in this paper. 



(20) 
(21) 
(22) 



2.4.1 NFW-EMD 



We convert the spherical NFW model to an EMD starting with the mass density defined by eqn. Q and letting r — >■ a„ 
rA — > HA, and 47r — >■ ATvp^qv, 



p(a) = 



Ma 



(23) 



-i-jvpyq^a^ A{c) ca (1 -I- caf ' 

where c = aA/ciscaic, and Oscaio is an ellipsoidal analog of r^. From this modified density we compute the mass enclosed within 
a„ by direct integration of p(a^) which gives. 



M(< a) = Ma 



$(x) 



A{c) 
le ex 
GMa c 



ln(l -f-ca) 



1 + ca 



Finally, from the expression for the density we construct the potential. 



du 



-GMa- 



ds 



2 A{c) A{u) 1 + ca{x, u) A{c) J,, S{s) 1 + ca{x, s) ' 

For Pv ~ qv ~ i this expression for $ reduces exactly to the spherical case (eqn. [3}. 



(24) 



(25) 
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2.4.2 CORELOG-EMD 

The procedure to construct the EMD analog of the spherical CORELOG model is fully analogous to that used for the 
NFW-EMD model. Converting eqn. ([SJ to an EMD gives the mass density, 



p{a) = 



Ma 



(1 + 



2. 3+(m)^ 



[1 + [caf] 



21 2 ' 



(26) 



where c = ai^ja^, and ac is an ellipsoidal analog of Tc. From this expression for the density follow the equations for the 
enclosed mass and the potential as a function of ellipsoidal radius, 



M(<a)^MA(l + c^) 

1 + (ca) 



and. 



GA/a 1 + 



du 



-GA/a 



1 + 



ds 



K^) 



iln 

2 



-In 

2 



1 + 



1 + \ca(x, 



l + [ca(f,u)]2 l + c2 

1 1 



1 + [ca(f, s)]2y l + [ca(f, s)]2 



1 + C2 



(27) 



(28) 



(29) 



Because the integral over a diverges at oo, we replaced the upper limit with a = 1. With this choice, we have for = g„ = 1 
that $ reduces to the spherical case (eqn. [7]) but offset by a constant so that, $(rA) = —GM/\/r^^, rather than "I>(rA) = 0. 



3 FIDUCIAL MASS MODEL PARAMETERS AND AXIAL RATIO PROFILES 



The ellipsoidal NEW and CORELOG models have four free parameters that need to be specified: concentration (c), mass 
(Ma), and axis ratios {jpv,qv). The axis ratios refer to the potential for an EP and to the mass distribution for an EMD. As we 
show below i n j}5l our results ar e not very sensitive to the degree of triaxiality of the ellipsoid as quantified by the triaxiality 
parameter T (jFranx et al.lll99ll ). Consequently, unless stated otherwise, throughout this paper we set. 



Pv 



(30) 



where we adopt T = 0.5 for the "maximally triaxial" ellipsoid that lies midway between the oblate (T = 0) and prolate 
(T = 1) spheroids. 

Since we also find that our results are very insensitive to the halo mass scale over the range of interest (10^^ — 10^'' Mq), 
we adopt a fiducial mass. Ma = 10" Mq, which lies near the middle of this range. With this value for the mass, we set the 
NEW concentration to c = 9 (assuming A ~ 10 0) appropriate for the results we obtained from X-ray measurements spanning 



the same large mass range IjBuote et al.ll2007^ . Although this normalization of the concentratio n-mass relation is 30-60% 

20081 ) . the results we present in 



Maccio et al 



larger than that obtained from dissipationless cosmological simulations (e.g.. 
this paper are unchanged if instead we use such a smaller concentration value. These values for c and Ma are used for both 
the EP and EMD variations of the NEW model. To provide a clear contrast with the NEW models, we adopt a large value 
of the concentration parameter (c = 100) for all CORELOG models which insures nearly scale-free behavior (e.g, M ~ for 
EP) over most of the system. 

Using these fiducial parameters for the EP models, we display in Figure [T] the (short-to-long) axial ratios qp of the 
isodensity contours plotted as a function of Uv/aA (with A — 100) for different values of g# = g^, the axial ratio of the 
potential e Uipsoid. We consider the r ange, g$ = 0.7— 1, corresponding to approximately, qp — 0.4— 1, appropriate for ACDM 
halos (e.g., Bailin fc Steinmetj 2005 . and references therein). The qp were computed for a given a„ by solving the equation, 
p(a„, 0, 0) = p(0, 0, z), for z and then setting qp = z/a^. We find that for both models, as is typical for EPs, that qp decreases 
(i.e., the density distribution becomes fiatter) with increasing a„. As approaches oa, qp approaches a constant value for 
each g$ for CORELOG-EP, whereas qp continues to decrease for NFW-EP; i.e., the NEW-EP mass distribution is fiatter 
than CORELOG-EP for a„ ~ oa and is slightly rounder for a„ ~ O.OIoa. (Note that for q,s> < 0.8 the density of the fiducial 
NEW-EP model takes negative values on the z axis for some z < oa-) We mention that the decrease of qp with radius 
for the NEW-EP model disagrees with the ge nerally increasing trend found in dissipationless cosmological simulations (e.g., 
Bailin fc Steinme tz!'2005': 'Allsood et al."2006"), although qp profiles that decrease with radius tend to result when dissipation 
is considered (e.g., Kazantzidis ot al. 2004; Dobattista et al. 2008). 

For EMD models with the fiducial parameters, we display in Figure [2] the axial ratios (g$) of the isopotential contours 
plotted as a function of a„/aA for different values of qp = q^, the axial ratio of the mass (density) ellipsoid. The q<s> were 
computed for a given by solving the equation, 4>(a„,0,0) = $(0,0,^), for z and then setting 5$ = z/a^. For both the 
NFW-EMD and CORELOG-EMD models 5.5 increases (i.e., the potential becomes rounder) with increasing radius, though 
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Figure 1. Axial ratio profiles of the mass density for the fiducial NFW-EP (Left PaneJ) and CORELOG-EP {Right Panel) models (ijaj. 
Shown are results for six values of the axial ratio in the potential . Note that for hydrostatic equilibrium the axial ratio of the X-ray 
emissivity is the same as the potential, independent of the temperature profile (see i|4.1l l. 
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Figure 2. Axial ratio profiles of the potential for the fiducial NFW-EMD (Left Panei) and CORELOG-EMD (Right PaneJ) models (@. 
Shown are results for six values of the axial ratio in the mass density qp. Note that for hydrostatic equilibrium the axial ratio of the 
X-ray emissivity is the same as the potential, independent of the temperature profile (see i|4.1ll . 



for CORELOG-EMD g<j attains a constant value for a„ ^O.laA for all qp explored. Generally, for a„ ~ oa the potential is 
rounder for NFW-EMD compared to a CORELOG-EMD model with the same qp. The models have very similar values of g$ 
for ~ O.OlaA- 
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4 METHOD TO TEST SPHERICAL AVERAGING 

Our objective is to quantify average biases and scatter in cluster properties arising from inconsistent spherical averaging 
procedures employed by observers and theorists, with particular attention devoted to properties measured by X-ray and SZ 
studies. It is most common for the theorist to spherically average the three-dimensional cluster available in their models, from 
which "true" spherically averaged quantities are derived. The observer, however, has access only to the cluster projected onto 
the sky. Generally, the observer circularly averages the cluster data projected onto the sky and then performs a deprojection 
assuming spherical symmetry. We refer to cluster properties inferred from such a deprojection as "observed" spherically 
averaged properties. By comparing the "observed" and "true" values for clusters of different shapes and projection orientations 
we quantify any bias due to spherical averaging in this way. Below we provide more details on how we make this comparison. 



4.1 Hydrostatic Equilibrium — Gas Mass and X-ray Emissivity 



To associate X-ray emission with a cluster entails filling the gravitational potential well with hot plasma, which is approximated 
very well by a non-rotating ideal gas in hydrostatic equilibrium, 

VPgas = -PgasV$, (31) 

where Pgas is the thermal pressure and Pgas is the gas density. For an isothermal gas, this equation can be solved for the gas 
density to yield, 

Pgas (5) = pgas.O exp (^^^ [$o " ^{x)] 



(32) 



where Pgas,o and "I>o are, respectively, the central values of the gas density and the potential, and T is the (constant) temperature 
of the gas. Our study focuses on the isothermal solution because temperature gradients have little impact on the flattening 
of X-ray images of galaxies and clusters. This is a consequence of the isopotential surfaces havin g the same shapes as the 
constant surfaces of gas density, temperature, and X-ray emissivity; i.e., the X-ray Shape Theorem ( Buote &: CanizaresI 1994 . 
1996al : iBuote fc Humphrevlboid ) . In fact, in Theorem 7 of Paper 1 we show that the observed deprojected spherically averaged 
mass profile for any EP model is independent of the temperature profile. Nevertheless, because other quantities may be more 
affected, we also explore models with temperature gradients using a p olytropic equa tion of state Pgas oc p^as- Inserting this 
relationship between the pressure and density into eqn. ((5T]) gives (e.g., Sarazin 198d ) ■ 



T{x) =To + 



7-1 ^ma 



s{x) = pgas.O 



7 fcs 
T{x) 



($0 -$(£)), 



To 



(33) 



(34) 



We focus our attention on models with 7 = 1.2 which best represent the fa lling temperature pr ofiles (outside of any cool core 



Ostriker et al.ll2005l ') 



of observed clusters and those produced in cosmological simulations (e.g., 

To define completely these models we need to specify the normalizations of the temperature and densi ty profiles of the 
hot gas. We normalize the temperature using the Afsoo — 2x relation given by Arnaud et al. 1 20051 . 2007 ). For isothermal 
models, we simply set T = Tx. For the polytropic models we determine the normalization constant To by computing the 
gas- mass- weighted temperature within a radius of vsoo and then set the result equal to Tx. The results presented below in ij5] 

are very insensitive to the adopted form of the M500 — Tx relation. 

Similarly, we normalize the gas density using the M500 — Mgas, 500 relation given bv lArnaud et al.l (|2007l ). That is, we 
determine Pgas.o by computing the total gas mass within a radius of rsoo and then set it equal to Mgas, 500 ■ While this 
normalization is convenient, it is not strictly necessary because we choose to neglect the self-gravity of the gas in our models 
("l-gas <'!')- see beginning of ^ 

Given the gas density and temperature we compute the X-ray e missivity, ex = p^ „.^MT), where A(r) is the plasma 



_t^a 

emiss ivity for a coronal plasma, which we take to be the APEC model (|Smith et al.ll200H ) as implemented in xspec (I Arnaud 



19961 ) integrated over energies 0.5-8 keV. We assume a constant metallicity of 0.5 solar throughout the system, though our 
analysis is not very sensitive to the assumed abundance profile. For isothermal models, we set A = 1 so that the emissivity is 
simply, ex 



2 

Pgas ■ 



4.2 True Spherically Averaged Profiles of Mass and other Quantities 

For a given mass model (^ we compute the true spherically averaged profile of the (total) mass as follows. For EP models the 
spherical averaging of the mass profile is conveniently achieved by applying Gauss's Law over a sphere of radius r (Theorem 2 
of Paper 1). There is no such simplification for the EMD models, where instead we directly integrate the density (i.e., either 
eqn. l23l or l26p over a spherical volume. We evaluate the spherically averaged mass profile in 500 (200) logarithmically spaced 
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bins for EP (EMD) models within a radius r = oa, where A — 100 (13]). By fitting this profile with the spherical version of the 
mass model (e.g., NFW - eqn.[2ll, we obtain the "true" spherically averaged parameters for the mass model, Csp and Ma,sp. 
For consistency with the observed mass analysis (see i]4.5p . the fitting is performed only within r = {pv^vy^^a^oo which is 
« rsoo. Note in most cases at least 80% of the logarithmically spaced bins are within this radius. Moreover, the inferred values 
of Csp and Ma,sp are not very sensitive to small changes (~ 10%) in the outer radius used in the fits. 

Having chosen both a mass model and a gas equation of state ( i]4.ip we compute the spherically averaged profiles of 
quantities directly associated with the hot gas. For EMD models we compute the "true" spherically averaged gas mass (A/gas, ap) 
by integrating pgas over a spherical volume of radius r = rsoo, where rsoo is evaluated using Afsoo, sp determined by fitting the 
mass profile as described above. For EPs the eflfective spherical average is used by setting Mgas.sp ~ Afgas(< {Pvqv)~^^^rsoo), 
which is the three-dimensional analog of Definition 1 of Paper 1. These procedures insure a consistent comparison of the gas 
mass with the observations ( >i4.5p . It then follows that the spherically averaged gas fraction is, /gas,sp ~ -A^gas,sp/Msoo, sp- For 
the emission weighted temperature (Tx.sp) we integrate exT within rsoo a nd then divide by the l uminosity within that volume. 
Using the temperature and gas mass we compute, lx,sp ~ A/gas, sp2x,sp ( Kravtsov et al. 20061 ). Finally, since the integrated 
Compton-y parameter is proportional to the volume integral of the gas pressure (e.g., eqn. 48 of Paper 1), we set Vsz,sp equal 
to the integral of Pgas within rsoo, ignoring proportionality constants which are unimportant for our study. 



4.3 Projection 

We construct the X-ray image by projecting the X-ray emissivity onto the sky for a given orientation. In Paper 1 we present 
analytical formulas for EPs relating the "observed" deprojected spherically averaged quantities (e.g., mass, gas mass, etc.) to 
their three-dimensional ellipsoidal counterparts. This means that for a given shape and projection orientation one may, e.g., 
immediately evaluate the "observed" deprojected spherically averaged mass profile for NFW-EP using eqn. (|12|) inserted into 
the result of Theorem 7 of Paper 1. Nevertheless, in our present investigation we find it useful to carry out the projection 
(and subsequent deprojection) numerically in order to test the accuracy of other aspects of our computer code, particularly 
via comparison to the EMDs which cannot be evaluated with simple analytical formulas. 

Since the X-ray emissivity of an EP is an ellipsoid, we follow the procedure described in § 5 of Paper 1 and carry out the 
projection by partitioning the emissivity into a series of concentric, similar triaxial ellipsoidal shells and by approximating 
the emissivity as a constant within each shell. For such a system viewed with orientation {8,(j>), expressions for the surface 
brightness Ex,i and temperature map {T}i are given by eqns. (66) and (69) in Paper 1. Each quantity is defined within 
an elliptical annulus i on the sky with inner semi-major axis a^-i, outer semi-major axis ai, and axial ratio qs (eqn. 18 of 
Paper 1) which is same for all annuli. 

We compute the surface brightness of an EMD, Ex(2;',2/') — J €x(x')dz' , via full integration of the X-ray emissivity 
along the line-of-sight for a given {9,<j)). Using the same procedure we compute the emission- weighted temperature map, 
{T){x',y') = J ex{x')T{x')dz' fExix' ,y'). The integrations proceed interior to the isopotential surface that passes through 
Ha on the major axis. Since the potential of an EMD is not itself ellipsoidal in shape, but is very nearly so, we approximate 
this boundary surface by an ellipsoid in our computations. 

We make no further modifications to the model surface brightness and temperature map. In particular, we do not fold 
the image through an image point spread function, nor do we fold the spectrum through a response matrix. We assume the 
observer is able to account for such issues perfectly, so our analysis is able to focus on potential biases arising only from 
geometrical issues associated with projection and spherical averaging of the galaxy or cluster. 

(Note that the sky map of the Compton-y parameter is obtained in the same manner as Ex by replacing the emissivity 
with the gas pressure.) 



4.4 Spherical Deprojection 

Now turning our perspective to that of an observer, we begin with the surface brightness, temperature, and Compton-y maps 
on the sky as produced in the previous section. In general these maps have non-circular shapes, but it is the goal of the 
fictitious observer to infer spherically averaged profiles from these maps. To this end we follow standard practice and bin 
the surface brightness into a series of A'' concentric, circular annuli defined so that Ri is the outer radius of annulus i. For 
consistency, we use the same bin definitions adopted in ^4.21 for the three-dimensional radius of the true mass profile; i.e., we 
associate each Ri on the sky with the equivalent r^ in three dimensions defined in M.2\ 

The procedure to create the circular binning diflers for the EP and EMD models. For EPs the constructed image is already 
binned in elliptical annuli, where the annuli definitions correspond to elliptical versions of the spherical radii adopted in ^4.21 
for the true mass profile; i.e., for each bin we have Oi = a„,,; — ri, where oii is the outer semi- major axis of elliptical annulus 
i on the sky and a„,i is the outer semi-major axis of the shell i in three dimensions. To achieve the circular binning on the 
sky, for Ex,i defined in each elliptical bin we associate an inner radius Ri-i = ai-i^/q^ and an outer radius Ri = a.iyjql (see 
Definition 1 of Paper 1). In contrast, for the circular binning of the EMDs we must resort to computing, J Ex (x', y')dx'dy' , the 
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full integration of the surface brightness over each circular bin {Ri-i, Ri), a much more computationally expensive procedure 
than for the EPs. To improve computational speed for the EMDs we use only ten points, equally spaced in azimuth, for each 
circular bin. We also take a dvantage of the el liptical s ymmetry and perfo rm the calculation over half of the circle. Finally, 
as in previous studies (e.g., Lewis et al. 20031 . see also McLaughlin 19991 ). we find it useful for each annulus to choose an 
intermediate average radius to represent the result for the entire bin. 

We obtain the observed radial volume emissivity in each shell i (e^ '"^") by depro.jecting the circularly a veraged surface 
brightness (S3f'J^) using the standard spherical "onion peeling" method ( Fabian et al. 1981 : Kriss et al. 19831 ). In particular, 
we convert the binned surface brightness values into luminosities, L"''^ = T^{Ri — Ri-i)^'x'i and then apply the onion peeling 
method as given by eqn. (66) of Paper 1 but using the spherical version of the projection matrix Vlj^^ . This procedure is also 
used for the SZ map where the pressure takes the place of the emissivity, and L'^"'^ is replaced by the Compton-y parameter 
within bin i. 

The observed spherically averaged gas density follows immediately from the emissivity, p^^^^"^^ ~ yf^^^, for isothermal 
models. For the models with temperature gradients, however, we require the temperature profile in order to evaluate the 
density. To deproject the temperature map we follow the same procedure as for the surface brightness with the following 
substitutions to eqn. (66) of Paper 1: Li — >■ {T}iLi and e — >■ eT. The result is the radial profile of the product of the 
emissivity and the temperature (^(T)^^'"^^ in three dimensions. Using this, along with the deprojected emissivity, gives the radial 



temperature profile, 



isp,obs 



(er)J^'°''°/e;^'°''°, which allows the gas density to be computed, p^J^°^'' 



(We assume the observer has perfect knowledge of the constant metallicity for evaluation of A. 



4.5 Observed Spherically Averaged Mass Profile 



From the observed spherically averaged density and temperature profiles, we compute the mass profile assuming hydrostatic 
equilibrium. 



M(< r)°p^= = 



rkBT° 



din p° 



dlnr 



+ 



dlnr° 



dlnr 



(35) 



where we evaluate the derivatives by interpolating the binned profiles with cubic splines (except for the outermost bin where 
we use logarithmic interpolation). We compute the mass profile in the same radial bins as done for the true mass profile ( i]4.2p . 
and then fit this observed mass profile with the same spherical version of the mass model to obtain the observed spherically 



averaged model parameters. 



and MX |p. In order to minimize any weak sensitivity of our results to the chosen outer 



boundary defined by the elliptical radius aioo in the potential (mass) for EPs (EMDs), we performed the fitting of the observed 
mass profile within a smaller radius. As remarked already in ^4.21 we adopted rsoo for this purpose, although we found little 
difference when using other conventional choices (e.g., r25oo, '"200)- The procedure to determine the other gas parameters is 
the same as described for the "true" values in ii4.2l For example, for both EP and EMD models we compute the spherically 
averaged gas mass (M°j 



by integrating /5gal,sp within a spherical volume of radius r = rsoo, where rsoo is evaluated using 



m; 



■obs 

500, spi 



from which immediately follows the spherically averaged gas fraction, /g. 



A^gas^p/Ai'soo, sp- In the same way we 



compute ls°z,sp within rsoo, though we note that in practice SZ observations without corresponding X-ray data must estimate 
rsoo in a different manner. 



5 RESULTS 

5.1 Preliminaries 

For our comparison of the "observed" ( ^4.5[) with the "true" ( ^4.2p spherically averaged quantities, we will focus on the 
percentage bias for each parameter; e.g., the total mass, 

where we suppress the notation for the spherical averaging, over-density, and "obs" denoting the observed quantity. With 
regard to the numerical accuracy of the biases presented below, there are differences between the EP and EMD models. In 
Paper 1 we showed that for EPs the deprojected spherical average of a parameter (e.g., MgJo, spl3 can be evaluated directly 
from the ellipsoidal distribution (e.g., M{< a„)) without the need to perform explicitly any projection and deprojection. 
Consequently, the biases for EPs can be evaluated extremely accurately and efficiently. Recall in i )4.3l we chose, nevertheless, 
to numerically project and deproject the EPs to assess the performance of our computer code. We found no significant 

^ For an EP this is also denoted by (M{< r^oo))'^'^ using the notation of Paper 1. 
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Table 1. Polynomial Approximations of the Average Biases as a Function of g„ 
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The orientation-angle-average bias and standard deviation as a function of qy approximated by 
a cubic polynomial: cq + ciQy + C2q^ + csq^. Results are given for the triaxial, isothermal NFW 
models evaluated at rsoo; c.f. Figures l4l [Sl Isl [Tol 

differences between the results obtained via numerical projection and deprojection with those obtained using the analytical 
formulas from Paper 1; i.e., the numerical uncertainty for the EPs is negligible for our investigation. 

However, the much higher computational expense of the EMD models, owing to the need for numerical integration of their 
non-analytical potentials, limits the number of radial and azimuthal grid points that can be practically employed, rendering 
the EMDs less accurate than the EPs. By experimenting with models having different numbers of radial and azimuthal grid 
points and having different tolerance precisions for the numerical integrations, we conclude that in most cases the numerical 
uncertainty of the EMDs is negligible for our investigation. The largest uncertainties usually occur for models with the smallest 
axial ratios (g„ < 0.5) which translate at most to bias errors of ~ 0.2% for EMDs. 

5.2 Bias Distributions 

We construct bias probability distributions by evaluating eqn. (|36p (appropriate for each parameter) for a large number of 
viewing orientations {6, (f>) generated by randomly sampling the solid angle over the orientation sphere. By taking advantage 
of the ellipsoidal symmetry of the cluster models, we restrict the orientation angles to the first octant of the sphere; i.e., 
< cos6 < 1 and < <jf> < tv/2. In Figure [3] we display the bias distributions generated by 10* random orientations of 
the isothermal NFW-EP model with g„ — 0.70 for the gravitational potential. All of the distributions are non-gaussian 
and somewhat asymmetric, though they are peaked toward the center. Interestingly, the mean bias for each parameter 
is significantly different from zero; i.e., biases are not eliminated through angle-averaging. As qv increases toward 1 the 
distributions of the isothermal NFW-EP model become progressively narrower and more symmetric (not shown). Despite the 
non-gaussian behavior for smaller q^, the mean and standard deviation provide a useful, convenient characterization of each 
distribution, which we employ below to frame our discussion of the systematic variations of the bias distributions with qv . 

5.3 Angle Averages and Standard Deviations 

In Figures IH O [H 1101 1121 and [13] we plot the mean and standard deviations of the bias distributions versus qv for the 
concentration, mass, gas mass, gas fraction, Ysz and Isz.x computed within rsoo for the isothermal models as well as the 
polytropic EMD models. The results for EP models reflect 500 random orientations while the more computationally expensive 
EMD models reflect 200 random trials. For comparison, for the case of the isothermal EPs we also display results for models 
computed within r25oo (Figure [14]) and for spheroidal NFW-EP models (Figure [15]) . Results for Tx and Yx using polytropic 
models are displayed in Figure [TT] All the figures plot the biases over the same range to facilitate visual comparison. Finally, 
in Table [1] we provide analytical approximations to the mean bias and standard deviation profiles for the triaxial isothermal 
NFW models computed within rsoo, which we find are represented very well by cubic polynomials. 

5. 3. 1 Concentration 

The concentration parameter displays the largest mean bias magnitude of the parameters investigated. While still not large in 
absolute terms, as seen in Figure [4] the average bias is nearly -2% for the smallest axial ratios considered. In addition, next to 
Ysz ( i|5.3.7p the concentration possesses the largest scatter. The NFW and CORELOG models have nearly identical average 
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Figure 3. Bias distributions generated from 10^ random trials in viewing orientation for the isothermal NFW-EP models with = 0.70 
for the gravitational potential. Displayed are results for the concentration (top, Jeft), mass (top, right), gas mass (middle, left), gas 
fraction [middle, right), Ysz {bottom, left), and Vsz.x (bottom, right). The vertical dashed lines show the mean of each distribution. 



biases, though CORELOG has a slightly larger scatter. Since the deprojected spherically averaged mass profile of an EP is 
independent of the temperature profile (Theorem 7 of Paper 1), temperature gradients do not affect the bias distribution 
for the concentration of an EP. Comparing the isothermal and polytropic EMD models in Figure 3] one sees that the EMDs 
are only a little affected by having different temperature profiles; i.e., the polytropic EMDs have slightly larger scatter and 
slightly smaller average bias magnitudes. The concentration biases are essentially the same when the results are computed 
within the smaller aperture r25oo (Fig. I14|l or when spheroids are adopted instead of the default maximally triaxial ellipsoid 
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Figure 4. Bias on tlie concentration parameter as a consequence of splierical averaging. (Left Panel) Results are shown for isothermal EP 
models plotted as a function of qv in the gravitational potential. (In hydrostatic equilibrium this is also the qv of the X-ray emissivity.) 
The solid lines are orientation angle-averaged values (NFW -white, CORELOG - blue), while the shaded regions (NEW - light blue/cyan, 
CORELOG - dark blue) are the la ranges. {Middle Panel) Same as the left panel except now the EMD models are shown and, therefore, 
qv corresponds to the axial ratio of the mass density profile. {Right Panel) Same as the middle panel except now results for the 7 = 1.2 
polytrope with a temperature gradient are shown. 



NFW 

CORELOG 



EP 



• Knv 

■ CORELOG 




EMD 




Figure 5. Same as Figure|4]except now the bias on the mass is shown. For clarity, the blue dashed lines also represent the la CORELOG 
region when it is contained within the la NFW region. Note that the la region for CORELOG-EP is negligible. 



(Fig, lisp . Perhaps the most notable difference is the slightly more negative average bias for the oblate NFW spheroids. Note 
also that because our implementation of the CORELOG model with a large concentration results in essentially a scale-free 
model, all the CORELOG parameter biases (including the concentration) are essentially unchanged when computed within 
''2500 compared to the default rsoo. 

Overall, the concentration parameter biases exhibit some general behavior that apply to most of the other parameters. 
These trends are usefully described with respect to the projections down each of the three principal axes, which we summarize 
here. Of the three principal axes, the smallest mean bias magnitude occurs for the intermediate-axis projection, while the 
projections along the short and long axes bracket the full range of biases for all orientations. The short-axis projection under- 
estimates the true value, and the long-axis projection overestimates it. For the specific case of the concentration parameter, 
these trends can be explained because the concentration is defined to be inversely proportional to the scale radius of the mass 
profile. It is reasonable to expect that the scale radius tends to be overestimated when the ellipsoid is projected along the 
short axis because the two longest axes then lie in the sky plane. Conversely, the scale radius tends to be underestimated 
when the ellipsoid is projected down the long axis because then the two shortest axes lie in the sky plane. 



5.3.2 Mass 

The NFW and CORELOG models exhibit interesting differences in their mass biases. The average bias of the NFW models 
is small but significant, increasing to ~ 1% for the smallest axial ratios (Fig. [Sjl. The la scatter for NFW models is similarly 
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0,6 

r/ r. 



Figure 6. Bias of the mass plotted as a function of radius between 0.25 — 1 rsoo for the isothermal NFW-EP model for = 0.7 {left 
panel) and = 0.9 (right panel) in the gravitational potential. Results are displayed for projections down the three principal axes: short 
(dashed, blue), long (dotted, red), and intermediate (solid, black). The top axis gives r/rs, where Ts is the NFW-EP scale radius. 



3 .1 2 3 4 




Figure 7. Bias of the mass plotted as a function of radius between 0.25 — Irsoo for the isothermal NFW-EMD model for q^ = 0.4 (Jeft 
panel) and g„ = 0.6 (right panel) in the mass. Results are displayed for projections down the three principal axes: short (dashed, blue), 
long (dotted, red), and intermediate (solid, black). The top axis gives r/r^, where is the NFW-EMD scale radius. 
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Figure 8. Same as Figure |4] except now the bias on the gas mass is shown. For clarity, the blue dashed lines also represent the la 
CORELOG region when it is contained within the la NFW region, and the average bias for CORELOG is denoted by the red solid line. 
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0,6 
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Figure 9. Bias of the gas mass plotted as a function of radius between 0.25 — Irsoo for the isothermal (solid lines) and polytropic 
(dashed lines) NFW-EP (left panel) and CORELOG-EP (right panel) models for qv = 0.7 in the potential. Results are displayed for 
projections down the three principal axes: short (blue), long (red), and intermediate (black). The top axis gives the radius expressed in 
units of the scale radius for NFW-EP and in units of the core radius for CORELOG-EP. 




small; < 1% for EP and about twice as large for EMD. Since the mass bias of an EP does not depend on the temperature 
profile, the results for the biases for the isothermal and polytropic EPs are identical (not shown). The strong similarity 
observed between the isothermal and polytropic NFW-EMD models in Fig. [5] indicates that the mass bias of EMDs is also 
mostly insensitive to the temperature profile. We obtain very similar results for the spheroidal NFW models (Fig. I15|) . 

In contrast, the isothermal CORELOG-EP model possesses essentially zero bias for all projections. This is a consequence 
of our adopting a large concentration (c = 100) so that CORELOG-EP approximates the scale-free logarithmic potential, 
which we demonstrated in Theorem 9 of Paper 1 (see also lChurazov et al.ll2008 ) always has zero mass bias when viewed from 
any orientation - independent of the temperature profile; i.e., the polytropic CORELOG-EP (not shown) also has zero mass 
bias for any projection. While not aways identically zero, CORELOG-EMD displays a very small range of biases (a < 1%), 
with the isothermal case having essentially zero average bias. The polytropic CORELOG-EMD model attains an average mass 
bias of about —0.4% for the smallest axial ratio studied (qv ~ 0.40), though we caution that the polytropic EMD models are 
the least accurate of all the models investigated (see end of ^5.1|) . 

For our fiducial models (excluding the zero-bias CORELOG-EP) the sign of the mass bias is opposite to that found for 
the concentration parameter ( i]5.3.1[) : i.e., projecting down the short axis results in a positive mass bias while the opposite 
happens when projecting down the long axis. We note, however, that if instead we adopt a smaller value for the concentration 
(say 10 instead of 100) for the CORELOG-EP model, the behavior reverses for the short and long axis projections. We 
conclude that the sign of the bias for the mass has a more complex origin tha n for the concen tration. 

Our results broadly agree with those obtained bv PifFaretti et al. ( 2003h and Gavazzil ( 20051 ) who have previously performed 
the most ext ensive ex a minat ions of the effects of spherical averaging on the the masses of galaxy clusters inferred from X-ray 
observations. Gavazzi (|2005l ) computed the mass bias as a function of radius for spheroidal, isothermal, NFW-EMD models 
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Figure 11. Similar to Figure |4] except now the biases on Tx and Yx are shown for the (Top Panels) polytropic EP models and {Bottom 
Panels) polytropic EMD models. Results for Yx for isothermal models are identical to those of the gas mass (Figure (8)1. For clarity, the 
blue dashed lines also represent the Icr CORELOG region when it is contained within the Icr NEW region, and the average bias for 
CORELOG-EMD is denoted by the red solid line. 




viewed down their symmetry axes. Visual inspection of Figure 6 of lGavazzil (|2005l 'l reveals that the mass bias for such face-on 
projections of spheroidal models is ge nerally only a few percent for a radius ~ 4rs ~ rsoo, which agrees reasonably well with 
our calculations. Piffarct ti et al.l ([2003) instead fitted ROSAT X-ray imaging data of ten clusters with an isothermal triaxial 
P model (i.e., isothermal CORELOG-EP) projected down the intermediate principal axis (i.e., edge-on). They obtained mass 
biases of magnitude < 4% broadly similar to those of our isothermal NFW-EP model, which provides a more appropriate 
comparison because our fiducial CORELOG-EP model has a much larger concentration than used by Piffarctti ct al. (2003|). 
When the NFW-EP model is computed within the smaller radius r25oo ~ 0.46r5oo (Fig. ll4p . both the mean bias (< 0.3%) 
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Figure 14. Results computed within ^2500 for isothermal EP models. 



and scatter (< 0.3%) are strikingly smaller than obtained for rsoo- Such large changes justify a more detailed examination of 
the radial dependence of the mass bias for the NFW models. (We do not include CORELOG in this discussion because for the 
large adopted concentration there are negligible changes in the mass bias with radius.) In Figure [6] we show the radial profile 
of the mass bias for NFW-EP models with = 0.7, 0.9. While the projection of the intermediate axis shows little radial 
variation, the other two projections have larger variations, each of similar magnitude, but of opposite direction. Interestingly, 
both the short- and long-axis projections cross the line of zero bias at nearly the same radius ~ 0.4 rsoo ~ l.Tr^. Although we 
have shown profiles only for two values of g„ , we find similar results for all g„ examined. 

For comparison, in Figure [7| we show the radial profiles for the mass bias of the NFW-EMD models with qv = 0.4,0.6, 
which have mass axial ratios at rsoo similar to the EPs shown in Figure [6] Although the intermediate-axis projection of 
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Figure 15. Comparison of results for oblate and prolate isothermal NFW-EP models. The solid lines are orientation angle-averaged 
values, while the shaded regions are the la ranges. Note the color coding is reversed for the gas fraction. 



the EMDs shows more radial variation than for the EPs, the characteristics of the profiles are very similar. In particular, 
the bias profiles of the short- and long-axis projections cross the zero bias line at nearly the same radius, this time near 
~ 0.36 rsoo ~ 1.5rs, which also applies to all the EMD models with different that we investigated. These results are little 
changed for the polytropic models, consistent with the insensitivity of the global mass bias to temperature gradients discussed 
above. The largest effects are observed for the = 0.4 EMD model where the crossing of the zero-bias line for the polytropic 
model occurs closer to 1.8rs, but the differences in the biases between the short- and long-axis projections at LSr^ are almost 
the same as at ~ I.Qts- 

Very similar results are also obtained for spheroidal models, where we have examined short- and long-axis projections, 
as well as projections along an intermediate direction defined by (9 — 45°, (j> — 45°). We conclude that a useful, general 
procedure to minimize the mass bias for a galaxy or cluster expected to follow the NFW profile is to measure the mass at a 
radius ~ 0.4r5oo ~ 1.7rs, which should be accurate to within « 0.5% for any intrinsic axial ratio (and triajdality) , projection 
orientation, and (isothermal or 7 = 1.2 polytropic) temperature profile. Since for our models r25oo ~ 0.46r5oo, the use of r25oo 
should also be a very suitable and convenient choice to minimize the orientation mass bias. 

As with the global mass bias, we find that our results for the radial variation of the mass bias broadly agree with the 
face -on projections of spheroidal isothermal NFW-EMD models examined by Gavazzi (2005). From visual inspection of Figure 
6 of iGavazzil (|2005l '). his symmetry-axis projections cros s the zero - bias l ine near a radius of 2rs, similar to the radius we find 
minimizes the mass bias for NFW models. Figure 6 of IGavazzil (|2005l ) also shows that the face-on prolate models tend to 
overestimate the mass within ~ 2rs and underestimate the mass for radii larger than ~ Sr^. The opposite is true for the 
oblate models. We confirm a similar trend for spheroidal models (not shown). 



5.3.3 Gas Mass 

The mean bias for the gas mass is very similar to that obtained for the total mass. As shown in Fig. [S] the largest bias is ~ 1% 
which occurs for the isothermal NFW-EMD model at q^j — 0.40, while all the CORELOG models have nearly zero average 
bias. (We obtain very similar results for the isothermal spheroidal NFW models shown in Fig. 1151 ) The bias scatter for the 
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gas mass is generally larger for the CORELOG models, and it exhibits a diversity of model-dependent behavior. For example, 
the presence of temperature gradients affects the scatter of NFW-EMD and CORELOG-EMD in the opposite manner; i.e., 
the scatter is larger for the polytropic model compared to the isothermal model for CORELOG and smaller for NFW. In 
contrast, for EPs the scatter is larger for the polytropic models for both the NFW and CORELOG (not shown). Finally, the 
sign of the bias in the gas mass for principal-axis projections is also model-dependent; e.g., for isothermal models CORELOG 
behaves as the concentration while NFW has the opposite behavior (like the total mass). 

When the isothermal EPs are computed within the smaller radius r25oo, the mean bias and scatter (Fig. I14p do not 
change by very much. The most notable difference is that the mean bias for NFW-EP is nearly zero for r25oo. It is also 
instructive to examine in more detail the radial dependence of the gas mass bias to compare to the total mass. In Figure [9] 
we display the radial profile of the gas mass bias for projections down the principal axes for the isothermal and polytropic 
NFW-EP models with q^j = 0.7. The bias in the gas mass for the isothermal model displays similarities to the bias in the 
mass, including crossings of the zero-bias line for the short-axis and long-axis projections near radius ~ 0.4r5oo. However, the 
polytropic model has qualitatively different behavior, showing larger biases that do not cross the zero-bias bias line. In fact, 
the polytropic NFW-EP profile is more similar to profiles of the CORELOG models (right panel of Figure [9| which display 
very little radial dependence as noted in H5.3.1I 



5.3.4 Gas Fraction 

The behavior of the mean and scatter in the bias of the gas fraction is extremely similar to the gas mass. Generally the 
magnitude of the mean bias of the gas fraction is the same as, or slightly smaller than, for the gas mass. 



5.3.5 Tx 

In Fig. [TT] we display results for the emission-weighted temperature integrated over the volume of default radius rsoo for 
polytropic models. Not only do the models all have nearly zero average bias, but overall the bias scatter is the smallest for all 
the parameters considered. We note that the sign of the bias (for principal-axis projections) for Tx behaves the same as for 
the mass. 



5.3.6 Yx 

For isothermal models the biases of Ix are identical to the gas mass since Fx = MgasTx. Given the extremely small bias and 
scatter for Tx ( ^5.3.5p . it is unsurprising that even for polytropic models the biases of Yx are extremely similar to the gas 
mass (Fig. Ilip . Hence, with regard only to the orientation bias of the type considered in our study, we find no significant 
advantage in terms of either the mean bias or scatter of using Yx instead of the gas mass alone. 



5.5.7 Ysz andYsz.K 



The integrated Compton-y parameter has a very small (^1%) mean bias for all the models considered (Figs. \T2l 1141 [T5)) . 
However, Isz does exhibit the largest bias scatter; e.g., for isothermal models, a = 6% for NFW-EP and a — 9% for 
CORELOG-EP. Although differing in detail, results for Fsz.x = Ysz/Yx, are fairly similar to those of Ysz. Probably the most 
noticeable difference is the smaller bias scatter for Ysz,x for CORELOG models. These results for Isz.x suggest that for small 
cluster samples attempts to calibrate the relationship between Ysz and Yx will encounter significant intrinsic scatter owing 
to orientation bias. Fortunately, our models indicate that by averaging over a large number of ran domly oriented clusters 



this source of systematic error c an be reduced to less than one percent. Current observations (e.g., lAndersson et al.l 12011 



Planck Collaboration et al.ll201ll ) are consistent with larger intrinsic scatter than that associated with the orientation bias on 
Ysz,x that we have computed. 



5.4 Comment on the Assumption of Random Orientations 

Our analysis assumes that clusters are randomly oriented to the observer's line of sight. However, in an X-ray flux-limited 
sample there may be a bias toward selecting clusters viewed down their major axis, since this orientation gives rise to the 
highest central surface brightness in an ellipsoidal cluster. As mentioned in fJT] we have assumed "ideal measurements" of 
galaxy clusters, which in this context means we have assumed a volume-limited cluster sample free of such observational 
biases. Similarly, the gravitational fields generated by large-scale structures likely establish weak cluster alignments on ~ 
40 Mpc scales (e.g., Paz ct al. 2011, and references therein). By assuming a sample volume much larger than this scale, we 
expect that small deviations from random orientations due to this effect impact our calculations less than the systematic 
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differences observed between the different models explored in our paper; e.g., EP vs. EMD, NFW vs. CORELOG, isothermal 
vs. polytropic. Cosmological N-body simulations are required to investigate this further. 



6 CONCLUSIONS 

This is the second of two papers investigating the spherical averaging of ellipsoidal galaxy clusters (and massive elliptical 
galaxies) in the context of X-ray and SZ studies. In Paper 1 we present analytical formulas describing both the deprojection 
and spherical averaging of EPs. In the present study we quantify the mean bias and scatter in cluster observables (e.g., mass) 
resulting from spherical averaging in the following sense. We associate "true" spherically averaged values with those computed 
from direct spherical averaging of the true three-dimensional cluster, as is done typically in theoretical studies. To obtain 
the observer's perspective, we initially fill up the ellipsoidal cluster's potential well with hot ICM by assuming hydrostatic 
equilibrium. The X-ray emission from this ellipsoidal model is projected onto the sky for a particular orientation. From this 
point onward the hypothetical observer treats the cluster as spherical; i.e., the observer circularly averages the X-ray image 
and temperature map and deprojects each of them assuming spherical symmetry. We use these profiles to compute "observed" 
spherically averaged quantities. By comparing the "true" and "observed" spherically averaged quantities we compute the bias 
arising from the inconsistent spherical aver aging procedures . 



Our study extends previous work ('e.g.. |Piffaretti et aklliooj : Icavazzilboosh on this topic in a few key respects. Of most 



importance, unlike previous studies that examined biases using only a single projection direction, we compute statistical 
orientation bias distributions for each cluster parameter. This allows us to present the first calculations of the orientation- 
averaged bias and scatter in cluster properties derived from X-ray and SZ studies. (Note that although we find that the 
parameter bias distributions are non-gaussian, because they are also centrally peaked it is nevertheless useful for our present 
investigation to describe them using the mean and standard deviation.) We also employ a more diverse set of cluster models 
than considered in previous individual studies. First, we consider clusters having a constant shape or one that changes with 
radius - in either the mass or ICM; i.e., both EPs and EMDs. Second, our fiducial models are maximally triaxial (triaxiality 
parameter, T = 0.5), though we also examine both oblate (T = 0) and prolate (T = 1) spheroids for comparison. Second, we 
examine gravitational potentials corresponding either to the NFW profile or a (nearly scale-free) logarithmic potential. Third, 
for the ICM we consider both isothermal models and (polytropic) models with temperature gradients. Finally, we examine 
the biases for several parameters: concentration, (total) mass, gas mass, gas fraction, emission- weighted temperature (Tx), 
Ix, integrated Compton-y parameter (i^sz), and the ratio Fsz.x = isz/Vx. 

We find that substantial biases can result from different viewing orientations, where Ysz and the concentration have the 
largest scatter (as large as a ~ 10% for Fsz) and Tx has the smallest (cr < 0.5%). As expected, the biases of largest magnitude 
occur for the fiattest models; i.e., = 0.7 for EPs and (jr„ = 0.4 for EMDs. Because the orientation bias of Tx displays such 
a small scatter, we find that the scatters for Ix and Afgas are virtually the same (as are their mean biases). Hence, spherical 
averaging of ellipsoidal clusters has essentially the same effect on Yx as it does on Mgas. We conclude in light of this that the 
smaller scatter for Fx obtained from cosmological hydrodynamical simulations ( Kravtsov et al. 20061 ) might be attributed to 



effects not considered in our study, especially deviations from ellipsoidal geometry and hydrostatic equilibrium. 

Although the mean biases are small and almost always < 1%, the ffattest models generally have a significant non-zero 
bias; i.e., typically orientation angle averaging does not completely eliminate projection biases. The mean biases of largest 
magnitude occur for the concentration parameter ~ —3%, while Tx always has nearly zero mean bias. Despite Ysz having the 
largest scatter, we find that its mean bias is always < 1%. Finally, the masses of the nearly scale-free CORELOG models - 
both EP and EMD - exhibit very small (or zero) mean bias and scatter. 

Orientation bias can generate significant scatter (up to cr ~ 8% for the ffattest systems) in the relationship between Ysz 
and Ix (as expressed by Vsz.x). This will be a significant source of uncertainty for measurements of the Ysz — Yx relation 
using small cluster samples. Fortunately, since the mean bias for Fsz,x is small (cr < 1%), for large enough cluster samples the 
orientation bias can be reduced below one percent. 

Our results do not change significantly over a wide range in halo mass (10^^ — 10^^ Mq) spanning massive elliptical galaxies 
to galaxy clusters. However, we do observe an interesting trend in the radial profile of the orientation bias for the mass for a 
cluster that follows the NFW profile. The scatter in the mass bias attains a minimum near the radius 0.4r5oo « l.lvs « r25oo- 
Consequently, considering the effect of orientation bias, the mass computed at this radius is always within « 0.5% of the 
true value for any orientation for all of our models. For comparisons between a small number of clusters, it is useful to quote 
masses at 7-2500 to minimize the orientation bias. 

It is important to recognize that even though the la scatter of the orientation bias for any parameter we have investigated 
is < 10%, the scatter is not small relative to other sources of systematic error. Take for example the gas fraction of the fossil 
group/cluster RXJ 115 9-1-5531 for which we pr ovide a detailed systematic error budget from a combined Chandra and Suzaku 



analysis (see Table 2 of lHumphrev et al.ll2011f ). The systematic errors reported on the gas fraction computed within rsoo are 



3% which is very similar to the la scatter in the orientation bias for the fiattest ellipsoids fFig. llOf) . Hence, orientation 
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bias is not negligible compared to other effects, and this is especially so in a fossil system like RXJ 1159+5531 where there is 
evidence that the hydrostatic equilibrium approximation is accurate. 

To facilitate the accounting for orientation bias in X-ray and SZ cluster studies, we provide cubic polynomial approxima- 
tions to the mean bias and la scatter as a function of axial ratio for each parameter for the (maximally triaxial) isothermal 
NFW models (Table [T|. These can be consulted to estimate the typical range of orientation bias for a given parameter and 
to construct a prior probability distribution for Bayesian stud ies. I n the latter case, to average over g„ for a given halo mass 

can be employed directly for the NFW-EMD model. 



the probability distribution for ACDM halos by Jing fc Sutol ( 



ACKNOWLEDGMENTS 

We gratefully acknowledge partial support from the National Aeronautics and Space Administration under Grant No. 
NNX10AD07G issued through the Astrophysics Data Analysis Program. 



REFERENCES 

Allen S. W., Evrard A. E., Mantz A. B., 2011, ARA&A, 49, 409 

AUgood B., Flores R. A., Primack J. R., Kravtsov A. V., Wechsler R. H., Faltenbacher A., Bullock J. S., 2006, MNRAS, 
367, 1781 

Andersson K., Benson B. A., Ade P. A. R., Aird K. A., Armstrong B., Bautz M., Bleem L. E., Brodwin M., Carlstrom J. E., 

Chang C. L., Crawford T. M., Crites A. T., et al. 2011, ApJ, 738, 48 
Arnaud K. A., 1996, in Jacoby G. H., Barnes J., eds, Astronomical Data Analysis Software and Systems V Vol. 101 of 

Astronomical Society of the Pacific Conference Series, XSPEC: The First Ten Years, p. 17 
Arnaud M., 2005, in F. Melchiorri & Y. Rephaeli ed.. Background Microwave Radiation and Intracluster Cosmology X-ray 

observations of clusters of galaxies, p. 77 
Arnaud M., Pointecouteau E., Pratt G. W., 2005, A&A, 441, 893 
Arnaud M., Pointecouteau E., Pratt G. W., 2007, A&A, 474, L37 
Baffin J., Steinmetz M., 2005, ApJ, 627, 647 
Binney J., 1981, MNRAS, 196, 455 
Binney J., Strimpel O., 1978, MNRAS, 185, 473 

Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition. Princeton University Press 

Borgani S., Finoguenov A., Kay S. T., Ponman T. J., Springel V., Tozzi P., Voff G. M., 2005, MNRAS, 361, 233 

Brighenti F., Mathews W. G., 2001, ApJ, 553, 103 

Buote D. A., Canizares C. R., 1994, ApJ, 427, 86 

Buote D. A., Canizares C. R., 1996a, ApJ, 457, 177 

Buote D. A., Canizares C. R., 1996b, ApJ, 457, 565 

Buote D. A., Gastaldello F., Humphrey P. J., Zappacosta L., Bullock J. S., Brighenti F., Mathews W. G., 2007, ApJ, 664, 
123 

Buote D. A., Humphrey P. J., 2011, MNRAS, in press (arXiv: 1109.6921) 

Buote D. A., Humphrey P. J., 2012, in D.-W. Kim & S. Pellegrini ed.. Astrophysics and Space Science Library Vol. 378 of 

Astrophysics and Space Science Library, Dark Matter in Elliptical Galaxies, p. 235 
Cavaliere A., Fusco-Femiano R., 1976, A&A, 49, 137 
Chandrasekhar S., 1987, Ellipsoidal figures of equilibrium 

Churazov E., Forman W., Vikhlinin A., Tremaine S., Gerhard O., Jones C, 2008, MNRAS, 388, 1062 

Debattista V. P., Moore B., Quinn T., Kazantzidis S., Maas R., Mayer L., Read J., Stadel J., 2008, ApJ, 681, 1076 

Evans N. W., 1993, MNRAS, 260, 191 

Evrard A. E., Metzler C. A., Navarro J. F., 1996, ApJ, 469, 494 
Fabian A. C, Hu E. M., Cowie L. L., Grindlay J., 1981, ApJ, 248, 47 
Fabricant D., Rybicki G., Gorenstein P., 1984, ApJ, 286, 186 
Franx M., Illingworth G., de Zeeuw T., 1991, ApJ, 383, 112 
Gavazzi R., 2005, A&A, 443, 793 

Grego L., Carlstrom J. E., Joy M. K., Reese E. D., Holder G. P., Patel S., Cooray A. R., Holzapfel W. L., 2000, ApJ, 539, 
39 

Henry J. P., 2003, in S. Bowyer & C.-Y. Hwang ed.. Astronomical Society of the Pacific Conference Series Vol. 301 of 
Astronomical Society of the Pacific Conference Series, Evolution of the X-Ray Properties of Clusters of Galaxies, p. 5 

Humphrey P. J., Buote D. A., Brighenti F., Flohic H. M. L. G., Gastaldello F., Mathews W. G., 2011, ArXiv e-prints 
(1106.3322) 



© 0000 RAS, MNRAS 000, 000-000 



Spherically Averaging Clusters: II. Biases 23 



Jing Y. P., Suto Y., 2002, ApJ, 574, 538 
Kassiola A., Kovner I., 1993, ApJ, 417, 450 

Kazantzidis S., Kravtsov A. V., Zentner A. R., AUgood B., Nagai D., Moore B., 2004, ApJ, 611, L73 

Krause E., Pierpaoli E., Dolag K., Borgani S., 2011, ArXiv e-prints 

Kravtsov A. V., Vikhlinin A., Nagai D., 2006, ApJ, 650, 128 

Kriss G. A., Cioffi D. F., Canizares C. R., 1983, ApJ, 272, 439 

Lewis A. D., Buote D. A., Stocke J. T., 2003, ApJ, 586, 135 

Maccio A. V., Button A. A., van den Bosch F. C, 2008, MNRAS, 391, 1940 

McCarthy I. G., Babul A., Bower R. G., Balogh M. L., 2008, MNRAS, 386, 1309 

McCarthy I. G., Schaye J., Ponman T. J., Bower R. G., Booth C. M., Dalla Vecchia C, Grain R. A., Springel V., Theuns 

T., Wiersma R. P. C, 2010, MNRAS, 406, 822 
McLaughlin D. E., 1999, AJ, 117, 2398 
Merritt D., Fridman T., 1996, ApJ, 460, 136 

Merritt D., Navarro J. F., Ludlow A., Jenkins A., 2005, ApJ, 624, L85 

Munoz-Cuartas J. C, Maccio A. V., Gottlober S., Button A. A., 2011, MNRAS, 411, 584 

Nagai D., Lau E. T., 2011, ApJ, 731, LIO 

Nagai D., Vikhlinin A., Kravtsov A. V., 2007, ApJ, 655, 98 

Navarro J. F., Frcnk C. S., White S. D. M., 1997, ApJ, 490, 493 

Navarro J. F., Hayashi E., Power C, Jenkins A. R., Frenk C. S., White S. D. M., Springel V., Stadel J., Quinn T. R., 2004, 

MNRAS, 349, 1039 
Ostrikcr J. P., Bode P., Babul A., 2005, ApJ, 634, 964 
Paz D. J., Sgro M. A., Merchan M., Padilla N., 2011, MNRAS, 414, 2029 
Piffaretti R., Jetzer P., Schindler S., 2003, A&A, 398, 41 
Piffarctti R., Valdarnini R., 2008, A&A, 491, 71 

Planck Collaboration Ade P. A. R., Aghanim N., Arnaud M., Ashdown M., Aumont J., Baccigalupi C, Balbi A., Banday 

A. J., Barreiro R. B., et al. 2011, A&A, 536, All 
Pratt G. W., Arnaud M., Piffaretti R., Bohringer H., Ponman T. J., Croston J. H., Voit G. M., Borgani S., Bower R. G., 

2010, A&A, 511, A85 

Rasia E., Ettori S., Moscardini L., Mazzotta P., Borgani S., Dolag K., Tormen G., Cheng L. M., Diaferio A., 2006, MNRAS, 
369, 2013 

Romeo A. D., Sommer-Larsen J., Portinari L., Antonuccio-Delogu V., 2006, MNRAS, 371, 548 
Sanders J. S., Fabian A. C, Smith R. K., 2011, MNRAS, 410, 1797 
Sarazin C. L., 1986, Reviews of Modern Physics, 58, 1 

Schuecker P., 2005, in S. Roser ed.. Reviews in Modern Astronomy Vol. 18 of Reviews in Modern Astronomy, New Cosmology 
with Clusters of Galaxies, pp 76-105 
Schuecker P., Finogucnov A., Miniati F., Bohringer H., Briel U. G., 2004, A&A, 426, 387 
Shaw L. D., Holder G. P., Bode P., 2008, ApJ, 686, 206 

Simionescu A., Allen S. W., Mantz A., Werner N., Takei Y., Morris R. G., Fabian A. C, Sanders J. S., Nulsen P. E. J., 

George M. R., Taylor G. B., 2011, Science, 331, 1576 
Smith R. K., Brickhouse N. S., Liedahl D. A., Raymond J. C., 2001, ApJ, 556, L91 

Takahashi T., 2011, in J.-U. Ness & M. Ehle ed., The X-ray Universe 2011, Presentations of the Conference held in Berlin, Ger- 
many, 27-30 June 2011. Available online at: http://xmm.esac.esa.int/external/xmm_science/workshops/2011symposium/, 
p.433 The ASTRO-H Mission, p. 433 

Tozzi P., 2007, in L. Papantonopoulos ed.. The Invisible Universe: Dark Matter and Dark Energy Vol. 720 of Lecture Notes 
in Physics, Berlin Springer Verlag, Cosrnological Parameters from Galaxy Clusters: An Introduction, p. 125 

Tozzi P., Norman C, 2001, ApJ, 546, 63 

Trinchieri G., Fabbiano G., Canizares C. R., 1986, ApJ, 310, 637 
Tsai J. C, Katz N., Bortschinger E., 1994, ApJ, 423, 553 

Urban O., Werner N., Simionescu A., Allen S. W., Bohringer H., 2011, MNRAS, 414, 2101 

Vazza F., Brunetti G., Gheller C, Brunino R., Briiggen M., 2011, A&A, 529, A17 

Voit G. M., 2005, Reviews of Modern Physics, 77, 207 

Voit G. M., Bryan G. L., Balogh M. L., Bower R. G., 2002, ApJ, 576, 601 

Voit G. M., Kay S. T., Bryan G. L., 2005, MNRAS, 364, 909 

White D. A., Fabian A. C, Allen S. W., Edge A. C, Crawford C. S., Johnstone R. M., Stewart G. C, Voges W., 1994, 

MNRAS, 269, 589 
White M., Hernquist L., Springel V., 2002, ApJ, 579, 16 
Younger J. D., Bryan G. L., 2007, ApJ, 666, 647 



© 0000 RAS, MNRAS 000, 000-000 



24 D. A. Buote and P. J. Humphrey 

Zaroubi S., Squires G., Hoffman Y., Silk J., 1998, ApJ, 500, L87 



© 0000 RAS, MNRAS 000, 000-000 



